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Abstract 

We study here the /i^/i^ decay mode of a very hght CP-odd Higgs boson of the 
NMSSM, ai, produced in association with a bottom-antibottom pair and find that, 
\ despite small event rates, a significant signal should be extractable from the SM back- 

■ ground at the LHC with high luminosity. 

cn ' 

1 Introduction 

In the Next-to-Minimal Supersymmetric Standard Model (NMSSM) |T], as a resuh of the 
introduction of an extra complex singlet scalar field, which only couples to the two MSSM- 
type Higgs doublets, the ensuing Higgs sector comprises a total of seven mass eigenstates: a 
^ . charged pair /i^, three CP-even Higgses /ii,2,3 ijnhi < < ""^/is) and two CP-odd Higgses 
^ ! ai,2 (^ai < ^02)- Consequently, Higgs phenomenology in the NMSSM is plausibly different 
from that of the MSSM, i.e, the minimal realisation of Supersymmetry (SUSY). 

In particular, over the past few years, there have been attempts to extend the so-called 
'No-lose theorem' of the MSSM - stating that at least one MSSM Higgs boson should be 
observed via the usual SM-like production and decay channels at the Large Hadron Collider 
(LHC) throughout the entire MSSM parameter space [2] - to the case of the NMSSM [21 HE]. 
From this perspective, it was realised that at least one NMSSM Higgs boson should remain 
observable at the LHC over the NMSSM parameter space that does not allow any Higgs-to- 
Higgs decay. However, when the only light non-singlet (and, therefore, potentially visible) 
CP-even Higgs boson, hi or /12, decays mainly to two very light CP-odd Higgs bosons, 
^1,2 — ^ aiOi, one may not have a Higgs signal of statistical significance at the LHC |6j. In 
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fact, further violations to the theorem may well occur if one enables Higgs-to-SUSY particle 
decays (e.g., into neutralino pairs, yielding invisible Higgs signals) [318]. 

While there is no conclusive evidence on whether a 'No-lose theorem' can be proved for 
the NMSSM at the LHC, there has also been put forward an orthogonal approach: to see 
when a, so to say, 'More-to-gain theorem' for the LHC [H [101 [HI 112] can be formulated 
within the NMSSM. That is, whether there exist regions of the NMSSM parameter space 
where more and/or different Higgs states of the NMSSM are visible at the LHC than those 
available within the MSSM. 

In our attempt to overview both such possibilities, we consider here the case of the yU+yU" 
decay channel of a light neutral CP-odd Higgs boson produced in association with 6-quark 
pairs at the LHC. This work complements the one carried out in which concentrated 
on ai —J- r+r" and 77 decays. 

The plan of this paper is as follows. In Sect. [2] we describe the parameter space scans 
performed. In Sect. [2] we discuss inclusive event rates. In Sect. [H we describe the signal- 
to-background analysis performed and introduce some benchmark points where an ai signal 
can be extracted via the decay mode. Finally, in Sect. [5l we summarise and conclude. 

2 Parameter Space Scan 

For our study of the NMSSM Higgs sector we have used NMSSMTools [131 [H]. This package 
computes the masses, couplings and decay widths into two particle final states of all the 
Higgs bosons of the NMSSM in terms of its input parameters, which can be taken at either 
the Electro- Weak (EW) or grand unification scale. NMSSMTools also takes into account 
theoretical consistency as well as experimental limits from negative Higgs searches at LEP 
[T5] and TevatronS, including the unconventional channels relevant for the NMSSM only. 

Here, instead of postulating unification or taking into account the SUSY breaking mech- 
anism, we fixed the soft SUSY breaking terms to a very high value, so that they give a 
small or no contribution at all to the outputs of the parameter scans. Consequently, we are 
left with six free parameters at the EW scale, uniquely defining the NMSSM Higgs sector 
at tree-level. Our parameter space is in particular defined through the Yukawa couplings A 
and K, the soft trilinear terms Ax and A^, plus tan/3 (the ratio of the Vacuum Expectation 
Values (VEVs) of the two Higgs doublets) and fies = A(S') (where (S) is the VEV of the 
Higgs singlet). 

In order to make a comprehensive study of the NMSSM parameter space, we have used 
the NMHDECAY code to scan over the aforementioned six parameters taken in the following 
intervals: 

A : 0.0001 - 0.7, k: 0-0.65, tan/3 : 1.6 - 54, 
/i : 100 - 1000 GeV, Ax : -1000 - +1000 GeV, A^ :-10 - 0. 

Soft terms which are fixed in the scan include: 
• = = rriD.^ = rriLs = ^E-i = 1 TeV, 

^Speculations of an excess at LEP which could be attributed to NMSSM Higgs bosons are found in [TO] . 
Very light CP-odd Higgs bosons of the NMSSM could also be produced in some rare hadron decays 
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. Au, = Ad, = Ae, = 1.2 TeV, 

• niQ = mu = rno = uil = uie = 1 TeV, 

• Ml = M2 = M3 = 1.5 TeV. 

In line with the assumptions made in [31 S], the allowed decay modes for the CP-odd 
neutral NMSSM Higgs boson ai are: 



We have performed our scan over 10 millions of randomly selected points in the specified 
parameter space. The output, as stated earlier, contains masses. Branching Ratios (BRs) 
and couplings of the NMSSM Higgses, for all the points which have passed the various ex- 
perimental and theoretical constraints. The points which violate the latter are automatically 
eliminated by NMSSMTools. 

3 Inclusive Event Rates 

The surviving data points are then used to determine the cross-sections for NMSSM Higgs 
hadro-production by using CalcHEP [ig and MadGraph [log, wherein some new modules 
have been implemented for this purpose [21] . As the SUSY mass scales have been arbitrarily 
set well above the EW one (see above), the NMSSM Higgs production modes exploitable 
in simulations at the LHC are those involving couplings to heavy ordinary matter only. 
Amongst the production channels onset by the latter, we focus here on 



i.e., Higgs production in association with a 6-quark pair. 

As an initial step towards the analysis of the data, we have computed the production 
cross-section times the decay BR against each of the six parameters of the NMSSM, as 
intimated. We started by computing total (i.e., fully inclusive) rates. Figs. 1 and 2 present 
the results of our scan, the first series of plots illustrating the distribution of event rates 
over the six independent NMSSM parameters, and as a function of the BR of the 
corresponding channel with the second plot displaying the correlations between the ai mass 
and the di-muon decay rate. It is clear that the large tan (3 and small fies (and, to some extent, 
also small A) region is the one most compatible with current theoretical and experimental 
constraints, while the distributions in k. Ax and are rather uniform (top six panes in 
Fig. 1). From a close look at Fig. 2 it is further clear that the BR(ai — )■ fi~^fi~) can be of 
0(10%) (0(1%)) [0(0.1%) or less] when 2m^ < rua, < 2m^(2m^ < < 2mb)[2mb < rua,], 
respectively. The first region of parameter space (m^^ < 2mr) is rather small, the second 
one (2mT- < nia^ < 2mi,) more significant, yet the widest one is the third one (2mb < ma J. 
However, by looking at the two bottom panes of Fig. 1, it is remarkable to notice that the 
event rates are sizable in all such mass regions, topping the 10'^ fb level in the two lower 

^We adopt herein CTEQ6L [20] as parton distribution functions, with scale Q — VS, the centre-of-mass 
energy at parton level, for all processes computed. 



ai — 7- yU^/i , r^r , gg, ss, cc, bb, tt, 77, Z7, sparticles. 



(1) 



(1(1,99 -^bb ai, 



(2) 
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mass intervals and the 10'^ fb level in the higher mass range. Finally, notice that the mass 
region below the //"""/i" threshold is severly constrained [22]. 

4 Signal-to-Background Analysis 

We perform here a partonic signal-to-background (S/B) analysis. We assume ^/s = 14 TeV 
throughout for the LHC energy. Also, in our numerical analyses, we have taken mb{mb) = 
4.214 GeV and m^"^'^ = 171.4 GeV for the (running) bottom- and (pole) top-quark mass, 
respectively, while we have input m^"^'^ = 1.777 GeV and m^"'*^ = 0.1057 GeV for the (pole) 
tau- and (pole) muon-lepton mass, respectively. After implementing the following standard 
cuts 

AR{b,b),AR{b,fx^),AR{b,fx^),AR{b,fi-),AR(b,fi-),AR{fx^,fx-) > 0.4 

|r/(6)|,|r/(6)|,|r/(/i+)|,|r7(/i-)|<2.5 

Prib), Prib) > 20 GeV, Pr(/i+), Prif^') > 5 GeV, (3) 

we obtain the invariant masses of the system depicted in Figs. 3-7, where we show 

the combined yield of the signal induced by qq, gg — )■ hbjX^ ii~ (via g and ai exchange) and 
of the irreducible background due to qq,gg — )■ bbfi^fi" (via g, 7 and Z exchange), including 
their interference. We also show in Figs. 3-7 the top-antitop reducible background, i.e., 
qq.gg ^ tt ^ bbW+W~ bbfi+ fi' P^'''. 

We notice the dominance of the 7 — )■ fj^~^fi~ tail of the irreducible background at very 
small di-muon invariant masses. Furthermore, the reducible background starts reaching its 
maximum at around Mw/'2. Finally, the Z — )■ fi~^fi~ peak of the irreducible background 
becomes overwhelming already at 60 GeV or so. Overall, the fj^~^fi~ signal yield in such a 
mass region is of order 75 (for reaching 60 GeV or so) to 1440 (for starting at 10 
GeV or so) signal events over a sizably smaller background (notice the logarithmic scale of 
the plots and notice that we are assuming, e.g., 300 fb~^ of accumulated luminosity). 

While the number of events is not very large, there is potential scope to extract a signif- 
icant signal over the above interval thanks to the high mass resolution that can be achieved 
using muon pairs. Assuming /i^/i" resolutions of 1 GeV [23j , we obtain as signal significances 
S/ \fB (where S and B are the signal and background rates, respectively, after a given lu- 
minosity), as a function of such a luminosity, those depicted in Fig. 8 (left hand side). The 
corresponding signal event rates are found instead in the right hand side of the same figure. 
From these last results, it is clear that detection at the LHC could occur for a\ masses be- 
tween 10 and ~ 40 GeV with rather modest luminosity, 30 fb^^ or so, while ~ 50 (~ 60) 
GeV masses require some 200(300) fb~^ while heavier states will not be resolvable even at 
the end of the LHC era. 

5 Conclusions 

Due to introducing a complex singlet superfield, the NMSSM can have a CP-odd Higgs boson 
with very low mass, a\. We have proven that there exist sizable regions of the NMSSM pa- 
rameter space where this kind of Higgs state, with a mixed singlet and doublet nature, could 
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potentially be detected at the LHC if 10 GeV < rua^ < 60 GeV in the ai — )■ mode if 

the CP-odd Higgs state is produced in association with a bb pair for large tan/3 plus small 
/ieff and A. After a realistic S/B analysis at parton level, we have in fact produced results 
showing that the extraction of light mass ai — fi'^fi resonances above both the irreducible 
and (dominant) reducible background should be feasible using standard reconstruction tech- 
niques 1211 [25]. While more refined analyses, incorporating parton shower, hadronisation 
and detector effects, are needed in order to delineate the true discovery potential of the LHC 
over the actual NMSSM parameter space, we are confident that our results are a step in the 
right direction to both: (i) prove the existence of a 'More-to-gain theorem' at the CERN 
collider for the NMSSM with respect to the MSSM (as fi~^fi~ signals from such light Higgs 
bosons are not at all possible in the latter scenario) and (ii) to establish a 'No- lose theorem' 
for the NMSSM at the LHC (as some of the parameter regions where the aforementioned 
signal can be detected overlap with those where hi^2 — ^ aiOi decays might be ineffective in 
extracting an hi signal). 

Notice that we have explored here the three mass regimes 2m ^ < rria^ < 2mr, 2mr < 
rUa^ < 2mb and 2mb < rua^ . The first interval in unresolvable because of a large irreducible 
background due to soft photons splitting into /x"*"//" pairs. The second one can actually 
be resolved only just below its upper end, i.e., close to the 2mi, threshold, above which ai 
state remains detectable for several tens of GeV. However, when m^^ > 60 GeV, it is the 
combination of both the irreducible and reducible background induced by on-shell Z bosons 
that prevents detection. 
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Figure 2: The CP-odd Higgs mass rria^ as a function of BR(ai — )■ ). 
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Figure 3: The differential cross section in the fj^~^fi~ channel for =9.76 GeV as a function 
of m^^, after the cuts in with A = 0.22341068, k = 0.4184933, tan/3 = 53.819484, /i = 
228.94259, Ax = -415.57365 and = -6.1773643. The solid line represents the signal and 
irreducible background together whereas the dashed line is the ti background. (Notice that 
here 2/71^°^*^ > nia^.) 
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Figure 4: The differential cross section in the channel for = 19.98 GeV as a function 

of m^^, after the cuts in © with A = 0.075946278, k = 0.11543578, tan/3 = 51.507125, fx 
= 377.4387, Ax = -579.63592 and A^ = -3.5282881. The sohd line represents the signal and 
irreducible background together whereas the dashed line is the tt background. 
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Figure 5: The differential cross section in the fi^fi~ channel for m^^ =30.67 GeV as a function 
of m^^, after the cuts in © with A = 0.10861169, k = 0.4654168, tan/3 = 48.063727, /x = 
222.99377, Ax = -952.59787 and = -7.2147327. The solid line represents the signal and 
irreducible background together whereas the dashed line is ti background. 
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Figure 6: The differential cross section in the jj'^jj channel for =46.35 GeV as a function 
of m^^, after the cuts in with A = 0.14088263, k = 0.25219468, tan/3 = 50.558484, = 
317.07532, Ax = -569.60665 and = -8.6099538. The sohd line represents the signal and 
irreducible background together whereas the dashed line is the ti background. 
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Figure 7: The differential cross section in the fi^fi channel for rria^ =60.51 GeV as a function 

0.17410656, K = 0.47848034, tan/3 = 52.385408, = 
= -9.0278415. The solid line represents the signal and 
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169.83139, Ay 
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455.85097 and A- 



irreducible background together whereas the dashed line is the tt background. 
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Figure 8: The significance S/ \fB (left) and total event rate S (right) of the qq, gg — ^ bbai — ^ 
bbfi^fi~ signal as a function of the integrated luminosity. 
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